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INTRODUCTION 


The  data  compiled  in  thi*  report  provides  engineering  information 
needed  in  making  evaluations  of  ammunition  improvements.  Represen- 
tative Caliber  . 50  and  20  mm  projectiles  data  appear  in  the  form  of 
computer  generated  plottings  of  trajectories  for  both  air  to  ground  and 
ground  to  ground  roles.  Specific  information  covers  time,  velocity 
and  altitude  versus  range.  Projectile  characteristics  include  muzzle 
velocity,  weight  and  shape  with  preformance  interests  centered  on  time 
of  flight,  muzzle  velocity  and  maximum  range  effectiveness.  Air  to 
ground  roles  include  fixed  wing  and  helicopter  launchings. 

'T  Caliber  . 50  projectiles  variations  listed  in 

Table  I comprised  a cross  section  of  nose  shapes  and  weights  for  the 
study.  Aside  from  potential  application  to  the  air  to  ground  and 
ground  to  ground  roles  these  projectiles  studies  have  accumulated 
quantities  of  experimental  information  needed  for  the  study. 


ITEM 

1 

2 

3 

4 

5 

6 

7 
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9 

10 

11 

12 

13 


Table  I 

Prolectlle  Studies 
WEIGHT 


PROJECTILE  61^S 


20mm  M56A3  1540 
20mm  M56A3  1950 
20mm  M56A3  2100 
20mm  PGU-3/B  1875 
20mm  PGU-3/B  1900 
20mm  PGU-3/B  2100 
20mm  M53  1540 
Cal.  .50  M8  662 
Cal.  .50  M8  850 
Cal.  .50  M8  1130 
Cal.  .50  M2  662 
Cal.  .50  M2  850 
Cal.  .50  M2  1130 


NOSE  SHAPE 

Blunt  - M505  Fuze 
Blunt  - M505  Fuze 
Blunt  . M505  Fuze 
Secant  Ogive  - SAPl 
Secant  Ogive  - SAPI 
Secant  Ogive  - SAPl 
Blunt  - M505  Fuze 
Tangent  Ogive 
Tangent  Ogive 
Tangent  Ogive 
Tangent  Ogive 
Tangent  Ogive 
Tangent  Ogive 
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METHOD 


Ir 
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TRAJE  (computer  program)  is  a point  mass  trajectory  program 
incorporating  variation  of  air  density  with  altitude  by  means  of  a sub- 
routine, Atmospheric  Computer  Program,  ATMSl,  (1).  The  basic 
TRAJE  model  (2)  is 


m 


dV 

dt 


(V)  - g 


m = Projectile  weight, 

V = Velocity  along  trajectory,  ft/sec 
p(y)  = Air  density  as  a function  of  height, 

A = Projectile  cross  section  area, 
g = Gravitational  constant  along  trajectory 
Cj^(V)  = Drag  coefficient  as  a function  of  velocity. 

Inputs  consists  of  projectile  weight  and  cross  sectional  area,  pro- 
jectile muzzle  velocity  and  drag  curve,  vehicle  velocity  and  angle  of 
ascent  or  descent.  Outputs  include  altitude,  range,  velocity,  time  of 
flight  and  drag  in  machine  plotted  form.  Output  information  is  used  to 
compare  trajectories,  time  of  flight  and  velocity  decay  of  various  pro- 
jectiles. 

Figure  1 contains  curves  (3)  of  drag  coefficient  versus  mach  num- 
ber for  projectile  items  number  1,4,7,  8,  and  11,  Appendix  A contains 
a listing  of  card  types  for  the  program  data  inputs,  and  the  units  for  the 
various  parameters.  The  program  is  written  in  Fortran  IV  for  use  on 
the  IBM  360  or  CDC  6500. 
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RESULTS 


A series  of  graphs  (Appendex  C)  present  the  charted  data  generated 
during  the  course  of  the  study  for  altitude,  time,  and  velocity  versus 
range  for  each  of  the  projectiles  investigated.  These  graphs  present 
various  combinations  of  the  same  basic  curves  to  provide  comparisons 
in  graphic  form.  Tables  II  to  V contain  listings  of  these  graphs.  Table 
II  covers  graphs  1 through  27  for  20  mm  projectiles  launched  in  the  air 
to  ground  role  at  a constant  aircraft  velocity  of  790  knots  and  a dive 
angle  of  30°  from  two  different  altitudes,  182.88  x 10  n^eters  (6000  ft) 
and  213.  36  X 10  meters  (7000  ft).  The  notations  in  the  column  headed 
"drag  Curve"  refers  to  the  M505A3  fuze  or  nose  shape  of  the  M56A3 
projectile  as  M505  and  to  the  nose  shape  of  the  20  mm  PGU-3/B  pro- 
jectile as  SAPI.  These  terms  also  appear  in  the  legend  of  the  plots  and 
in  Table  I and  are  used  to  identify  the  drag  curve  used  in  the  calculations. 

Graphs  1 and  2 compare  projectiles  of  two  different  weights  having 
the  SAPI  nose  launch  at  two  diferent  altitudes.  Graphs  3 and  4 com- 
pare projectiles  having  the  M505  nose  shape,  same  weights,  but 
different  muzzle  velocities  launched  at  two  different  altitudes.  Graphs 
5 and  6 provide  the  same  comparison  for  projectiles  having  the  SAPI 
nose  shape.  Graphs  7 through  12  compare  projectiles  having  M505  and 
SAPI  nose  shapes  at  two  weights  and  three  velocity  levels  launched  at 
two  altitudes.  Graphs  13  through  16  present  a comparison  of  the  M505 
and  SAPI  nose  shapes  at  two  launch  altitudes  by  interchanging  projectile 
weights.  Graphs  17  and  18  compare  SAPI  nose  shapes  at  two  launch 
altitudes  and  projectile  weights  and  three  velocities.  Graphs  19  and  20 
cover  the  same  comparison  for  the  M505  nose  shape.  Graphs  21  and 
22  compare  a projectile  having  the  standard  M56A3  weight  with  M505 
and  SAPI  projectiles  having  the  same  weights.  Graphs  23  and  24  pre- 
sent a comparison  of  five  projectiles,  SAPI  and  M505  projectile  nose 
shapes  at  two  different  weights  and  the  standard  M56A3.  Graph  25  com- 
pares the  standard  M56A3  weight  projectile  to  two  heavier  M505  and 
SAPI  nose  shape  projectiles.  Graph  26  compares  equal  weight  and  vel- 
ocity M505  and  SAPI  projectile  shapes  with  standard  M56A3.  Graph  27 
includes  the  M53  API  data  for  comparison  with  the  M56A,  M505  and 

SAPI  nose  shaped  and  weight  projectiles. 

Table  III  contains  a listing  of  Graphs  28  through  76  together  with 
identification  of  the  conditions  for  the  graph  data.  These  conditions 
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differed  from  those  identified  in  Table  II  by  varying  altitude  and  air- 
craft speed  at  lower  levels  and  dive  or  descent  angles.  The  range 
selected  relates  to  helicopter  launchings.  Four  projectiles  were  con- 
sidered; standard  M33,  one  with  M505  nose  shape,  and  two  SAPI  nose 
shaped  projectiles  differing  in  weight  and  velocity.  Each  graph  con- 
tains four  sets  of  data,  one  set  for  each  projectile.  Table  III  groups 
the  graphs  according  to  aircraft  velocity:  Graphs  28  through  37  are 
for  0 forward  velocity  while  Graphs  68  through  76  are  for  200  knots 
air  speed  with  the  other  graphs,  38  to  67,  grouped  for  velocities  in 
between. 

Table  IV  lists  24  graphs  covering  trajectory  data  for  20  mm  pro- 
jectiles in  the  ground  to  ground  role  at  elevations  in  the  3 to  60°  range. 
Projectiles  having  M505  and  SAPI  nose  shapes  at  two  velocity  levels  are 
compared  at  one  weight  level  together  with  the  standard  M56A1  weight 
projectile.  Graphs  77  through  89  cover  one  velocity  level  and  Graphs 
90  through  100  another  for  the  M505  and  SAPI  nose  projectiles. 

Table  V contains  a listing  of  the  graphs  presenting  Caliber  . 50  data. 
The  data  involves  two  Caliber  . 50  type  nosed  projectiles,  M2  and  M8^ 
launched  at  three  different  weights  over  an  elevation  angle  range  of  1 
to  60°.  Graphs  101  through  120  relate  to  the  M2  data  and  Graphs  121 
through  140  to  the  M8  data. 


Table  II. 

20»1  Projectllea,  Air  to  Ground  Trejectorlee 
Aircraft  Velocity  790  Knota.  Dive  Angle  30* 


Altleuda 

rt/Wet»re 


Kuecle  Telecltp 

tiilms 


1 

MOO/IS2S.I 

6AP1 

1950 

29S0/«».1« 

SAPI 

1850 

30U/919.2I 

2 

7000/2133.6 

SAPl 

1950 

2990 

1890 

9018 

3 

7000 

M505 

1990 

H505 

1990 

2950 

A 

6000 

M505 

1990 

2950 

M505 

1990 

9050 

5 

7000 

flAPI 

1990 

5050 

SAPI 

1990 

2990 

t 

6000 

SAPI 

1990 

2950 

&API 

1990 

5050 

7 

7000 

M50J 

1990 

2990 

SAPI 

1950 

2990 

t 

6000 

SAPI 

1990 

2990 

M509 

1950 

2990 

« 

6000 

SAPI 

1990 

5090 

M505 

1990 

5050 

10 

7000 

M505 

1990 

9050 

SAPI 

1950 

5090 

11 

7000 

M505 

aoo 

2B9O 

SAPI 

aoo 

2850 

12 

6000 

SAPI 

2100 

2890 

H905 

2100 

2890 

U 

7000 

SAPI 

2100 

2890 

M909 

1990 

2990 

M505 

1990 

3050 

lA 

6000 

SAPI 

2100 

2850 

M505 

1990 

5090 

H505 

1990 

2990 

13 

7000 

H909 

2100 

2890 

SAPI 

1950 

>050 

SAPI 

1990 

2950 

1« 

6000 

M909 

2100 

2890 

SAPI 

1990 

2950 

SAPI 

1990 

>050 

17 

7000 

SAPI 

2100 

2890 

SAPI 

1990 

5090 

SAPI 

1990 

2950 

IS 

6000 

SAPI 

2100 

2690 

SAPI 

1990 

>090 

SAPI 

1990 

2990 

IS 

6000 

H905 

aoo 

2850 

M909 

1990 

2990 

H909 

1990 

9090 

_ I 


i 

f 
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Table  II. 

- 

continued 

Gr^fh 

Nmn^T 

Altltu4« 

rc/N*MT« 

Dran  Curve 

Mltht 

Cnlna 

Mu«al«  Valeclty 

?0 

7000 

M505 

2100 

2850 

M505 

1950 

J050 

M505 

1950 

2950 

21 

7000 

«5°5(M56) 

15'*0 

3350/1021,08 

«505 

2100 

2850 

SAPI 

2100 

2850 

22 

6000 

«505(M56) 

1540 

JJ50 

M505 

2100 

2850 

SAPI 

2100 

2850 

23 

7000 

SAPI 

1950 

2950 

M505 

1950 

2950 

SAPI 

2100 

2850 

■505 

2100 

2850 

H505(M^6) 

1540 

3550 

2A 

6000 

M505(m56) 

1540 

3350 

M505 

2100 

2850 

SAPI 

2100 

2850 

«505 

1950 

2950 

SAPI 

1950 

2950 

23 

6000 

M505(M56) 

1540 

3350 

«505 

2100 

2850 

SAPI 

2100 

2850 

6000 

M505 

1950 

3050 

SAPI 

1950 

30J0 

26 

7000 

SAPI 

1950 

3050 

«505 

1950 

3050 

SAPI 

2100 

2850 

M505 

2100 

2850 

M505(M56) 

15‘tO 

3350 

27 

6000 

«505(N56) 

1540 

3350 

API 

2500 

2650/807.72 

SAPI 

2100 

2850 

«505 

1875 

3016 

SAPI 

1875 

3016 

♦NOTE  I: 

The  M505  Dr«g 

curve  was  obtained  froo  Ballistic 

Research 

Laboratorji  Aberdeen,  Hd 

and  SAPI  drag  curve  was 

obtained 

from  20mm  API 

and  SAPI  Design  Review  HcsoranduA,  Phsss  October 

2«,  1969, 

by  AVCO  Corporation 

Ordnance 

Dlvlalon. 

I 


f 


/ 
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Drag  Curve 
used 

TABLE  III. 

JOm  PKOJICTILES 
ADL  TO  GBOUND  TRAJICTORIES 

Weight 

grains 

Muzzle  Velocity 
f ps/raps 

M53 

1540 

3100/944.88 

SAPI 

2100 

2450/746.36 

SAPI 

2500 

2150/655.32 

M505 

2100 

2450/746.36 

Graph 

Altitude 

Helicopter 

Descent 

Number 

feet/meters 

Velocity  knots/fps/mps 

Ancle 

28 

200/60.96 

0.0 

5" 

29 

1000/304,8 

0.0 

5“ 

30 

200 

0.0 

10® 

31 

1000 

0.0 

10® 

32 

200 

0.0 

150 

33 

1000 

0.0 

15® 

34 

200 

0.0 

30® 

35 

1000 

0.0 

0 

0 

36 

200 

0.0 

45® 

37 

1000 

0.0 

45® 

38 

200 

50.0/84.45/57.58 

5* 

39 

1000 

50.0 

5“ 

40 

200 

50.0 

10® 

41 

1000 

50.0 

10® 

42 

200 

50.0 

15“ 

43 

1000 

50.0 

15“ 

44 

200 

50.0 

30® 

45 

1000 

50.0 

30® 

46 

200 

50.0 

45® 

47 

1000 

50.0 

45® 

48 

200 

100.0/168.89/115.  16 

5“ 

49 

1000 

100.0 

5“ 

50 

200 

100.0 

10® 

51 

1000 

100.0 

10® 

52 

200 

100.0 

15“ 

53 

1000 

100.0 

15“ 
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TABLE  III 


continued 


Greph 

Altitude 

Helicopter 

Descent 

Nue^cr 

feet/aetere 

Anele 

54 

200 

100.0 

30^ 

55 

1000 

100.0 

30* 

56 

200 

100.0 

45* 

57 

1000 

100.0 

45* 

58 

200 

150.0/253.34/172.73 

5* 

59 

1000 

150.0 

5® 

60 

200 

150.0 

10* 

61 

1000 

150.0 

10* 

62 

200 

150.0 

15® 

63 

1000 

150.0 

15® 

64 

200 

150.0 

o 

O 

65 

1000 

150.0 

30® 

66 

200 

150.0 

45* 

67 

1000 

150.0 

45* 

68 

200 

200.0/337.79/230.31 

5* 

69 

1000 

200.0 

5“ 

70 

200 

200.0 

10® 

71 

1000 

200.0 

10® 

72 

200 

200.0 

15® 

73 

1000 

200.0 

15® 

74 

200 

200.0 

30® 

75 

1000 

200.0 

30“ 

76 

1000 

200.0 

45® 

10 


TABLE  IV. 

20  MM  PORJECTILES 
^jHQUNDJTQ  ground  TRAJECTORTirf;* 


Graph 

Number 

Drag 

Curve 

Weight 

Grains 

Muzzle  Velocity 
fps/moB 

77  - 88 

M505  (M56) 

1540 

3350/1021. 08 

M505 

1950 

2950/899. 16 

SAPI 

1950 

2950/899. 16 

89  - 100 

M505  (M56) 

1540 

3350/1021. 08 

M505 

1950 

3050/929. 64 

SAPI 

1950 

3050/929. 64 

♦ 

ANGLES  OF 

ELEVATION:  3°, 

6°.  9°.  12°,  15°, 

18°,  21°,  24°,  27° 

30°. 

45°.  60° 

TABLE  V, 

CALIBER  . 50  PROJECTILES 

GROUND  TO  GROUND  TRATECTnnirc* 

Graph 

Number 

Drag 

Curve 

Weight 

Grains 

Muzzle  Velocity 
fps/mps 

101  - 120 

M2 

1033. 

2488/758. 34 

M2 

850. 

2910/886. 97 

M2 

662. 

2600/792. 48 

121  - 140 

M8 

662. 

2910/886. 97 

M8 

850. 

2600/792.48 

M8 

1130, 

2488/758. 34 

^ANGLES 

OF  ELEVATION; 

. o O O O 

1 , 2 , 2.  5 , 3 thru 

.,0^0  o o 

16  , 30  , 45  , 60 

11 
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CONCLUSIONS 


This  report  covers  a wide  range  of  combinations  of  drag  curve, 
projectile  weight,  and  muzzle  velocities.  These  parameters  in  turn 
effect  the  time  of  flight  and  maximum  range  of  the  performance  of  the 
various  projectiles.  It  is  recommended  that  similar  trajectories  be 
generated  for  an  air  to  air  role  as  well  as  other  ammunition. 
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TRAJE  DATA  INPUT  AND  FLOW  CHART 


CARO  TYPE 


COLUMNS 


FORMAT 


SYMfOL 


identification 


UNITS 


Read  catd  3 
If  carvi  8 E 


Calculates  tra'cctorv 
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APPENDIX  B 

nOMPUTER  program  LISTING 


PPOOHAM  TP*Jt( INPUT. OUTPUT. TAPE  I » I NPUT . T APE3*nuTPUT . T APE*. ) 

PEAL  LOT 

niMENSlON  SM«PE(h» 

01mEN'^10NANOA(200)  .CO (200)  .VCLUa(?00) 

COMMON  Y(?00).T(?00) .V(200) .CO (POO) .K.  I.X(POO) 

COMMON  VV(?00I .CC(POO) .VPLANE 
COMMON  27(?00) 

^ FOwma)  (6Ah) 

FOWMaT  (POX.  IPH  MU2ZLE  VELOCITY  « . FS.O.  7m  FT/SEC  ) 

7 FOPMAr  ( )M  . pA7) 

P format  (PFS.Pi 
P EOPmat  (PASI 

PS  FOPMAT  (IHl.  lOX.  70H  aNGO  (M)LS).  0(LBS/FT3).  PPES  (IN).  TEMP  (F) 
A XM  (fi-JAlNS)  VKET/SEC)  ) 

P7  FnKMAT(  PSh  «(1N0  velocity  (M/SEC)  « . F7.A.  1?H  PPtSbUPE  * .F7.<.. 
PlSM  TE'-PEH-iT  JMF(F)  * . F7.S) 

AA  FOP-IATdOX.  1 IH  Y(METt»S).  SiX.llM  X (METERS)  . 13X.  /H  T(SEC).  6X.  1 
Pah  V («ETEPS/SEC) . 16X.  ah  CO  . 9X.  11h  Z(HETERS)  ) 
aS  FOPMATdX.lFPD.rt) 

Ah  FOhMATdX.lOH  altitude  .lOX.  9H  DENSITY  . llX.  I 3m  TEMPERATURE  .7 
?X 

ss  fOwMAT(  SFPO.h) 

SS  FOh'MA  r (SFlfe.fl) 

S7  FOHMATdOOX.iFlO.A) 

SA  rOPMAT  dX.6F20.8) 

AO  FOkMATOIH  OEhSITY  » .F7.6.  8H  TEMP  r .FIO.A.IOH  HEIGHT  * .FIO.A) 
A)  E'OHMAT  (iHl) 

AE  FOhMAT(SFlO.S) 

AS  FORMAT)  isx.  aEIS.B) 

7A  FOkMAK  1AM  FORCE  (NTS)  e.ElA.6.  1 5H  YDOT  (H/SEC)  *.E16.6.  1 6M  VEL 
2 (M/SFC)  « .E16.b) 

EO  FORMATdAH  VEL(M/SEU  * . lAX  . 7H  CO  « . 20X) 
fll  FOk  (at ( I7.AF1S.8) 

10?  Frfh.MAT  dSFS.O) 

103  FOm mat  ( IHl  . Al A 1 1 
KlA  FOR MAT (F7.2) 

202  FOE  lAT (2F20.E) 

203  FORMaT(FA,2) 

20A  fomhat (Fa. 3. fa. 51 

POA  FOkMAT ( 3F7.A) 

A?«  format  (IHh  plane  velocity  ».  F6. 0.3AM  FT/  SEC  AT  AN  ANGLE  OF  OESC 
lENT  OF  . F6.A.  9H  RADIANS.) 

S?E  FORMAT (A6H  INPUT  IS  REL  PROJ  VEL  IN  FT/SEC. DRAG  COEF  KD.) 

1000  FOR  (AT  ( 3A 1) 

1010  FOrhaT(  ilH  wl  OF  PROJ.FIA.H.  IIH  NOSE  SHAPE. AS) 

1011  FormaT(13h  mij/zle  VEL*  . E1A.7.  ftH  FT/SEC  ) 
lOlP  E ORmaT (1 X.hAA ) 

ion  FORMATdX.ElA.fl.lRM  WT.  OF  PROJECTILE  ) 

’****•*••*  DATA  input  •******•••*•*****• 

1 PRtSSURE. temperature  (PF7.A) 

? SHAPF  (6AA) 

••••••«  3 VV(J) (FT/SEC). CC(J) (KO)  Fb.3.F8.5) 

A OVV( INCREMENTS/VELOCITY)  (F7.P) 

S Yd)  (AlTITUOE)  (F7.P) 

A G(GRAVITy-METERS  per  seconds  squared  .A(CROSSSECTIONAL  AREA) 
.XM(MASS-GRAInS) . ANGO (ANGLE  Of  DESCENT  OH  ASCENT . V 1 ( MUZZLE 
7 WIND  VELOCl  rr  (F7.A) 

M AIRCRAFT  VELOCITY  (Fb.2) 

P AiTlUDE  **  FOR  same  aircraft  VELOCITIES  CAN  RUN  AT  VARIOUS  ALTI- 
TUDES. TO  change  aircraft  velocity  Put  -lo.  follox 
XITM  altituue  and  then  aircraft  velocity  will  return 

ASKING  FOR  ALTITUOE  AGAIN  *• 

10  -10.  PLACED  FOR  LAST  AIRCRAFT  VELOCITY.  -10.  PLACED  FOR  LAST  G 
N*NUMHER  OF  TIME  INTERVALS 
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C OI*TlMt  INCKtMENTS  <StC» 

C VFLOr*DEC»E*SE  IN  VELOCITY  DUE  TO  0«*0/0»*G  VEL  IN  Y DIRECTION  (M/SEC) 
C »V**VER*GE  0ECELEM4T10N  IVEU  DUE  TO  DRAG  <M/SEC» 

C**************  ATHSI  « FOOIO 

C ATMSl <YH*ALTITUOE<rEET> . TEMPIRANKINE) .PRES  <SLUG/FT»*2. 

C DENSITY*  SLUG/rT«»J) 

C A*  CROSS-SECTIONAL  AREA  FOR  20RM  lEITMER  .515  OR  PI»R«R» 

C C*CONSTANT 

C ANG*  angle  of  TRAJECTORY  . YOOT*  Y VELOCITY  »OA*CHANGE  IN  X OISTANCE 
C X(I)*X  OISTANCE  • AV  * AVERAGE  VELOCITY  . V(I)*  Y OISTANCE 
C VEL*  FT/SEC  .KO  DRAG  CURVE 

C 7 IS  RiNO  VELOCITY  effect  ON  PROJECTILE 

c«*  temperature  in  F 
c**  pressure  C“  of  mercury 

c***«  V**  MINO  VELOCITY  IN  METERS  PER  SECOND 

C»*  IF  VFLOCITYIFP-,)  ANO  KU  CURVE 
C**  VVU)  *VV  < J) /3.2M0R 

201  CC(J)»CC(J>*H./3.IAI59 
IF  INPUT  VELOCITY (MACH)  ANO  CO  CURVE 


C*» 

C*« 

C* 

C«* 

C*» 


SPSO*SO'-T(<'.5R.G7  » TEMP) /5IB.6T)  •3A0.2RA 

VV(J)*VV(JI*SPSO 
201  CONTINUE 
IN*  I 
10*  T 

? CONT INUE 

WEAUC  IN.20#))P,TEHP 
REAO< IN.S) ShAPE 

SPSO*SORT  ( <<i5V.67  ♦ TEMP ) /5 16.67)  *3A0 . 2RA 

00  201  J*1.200 

RE»0<IN,204)VV<J).CC<J) 

IF(VV<J) .LF.-  9.)G0T022H 
1F(VV<  J)  .GF.IOOOO.)GUTU.tO 
VV(J) *VV<J)«SRSO 

201  continue 

M.J-I 
OV*VV  ( I ) 

N»VV  <1 ) /lO. 

F«I 

no2nj*2,N 

OV«OV-IO. 

call  1NTERR<CC.C0< J) .OV.VV.IO 
V( J) *ov 

20  CONTINUE 

V(  U «VV  < I )• 

CD<n»cc<n 

WRITEU  ).?5) 

KKK»0 

nn  1001  j*i.N 

VV< J) *V  <J) 

CC< J) »C0< J) 

1001  CONTINUE 

READ*  IN.  I06)0W 

21  CONTINUE 
REAUUN.IOM  YU) 

IF<YU).LE.-IO.)GOTOl7 

r.OTOlM 

17  CONT INUE 

REAnUN.lOG)  YU) 

REAOUN.203)  yPLANE 
IF«VPLANE.LE.-I0.)GUTO19 
Ifi  CONTINUE 

IF (KKK.GF. 1 )GUT02m 
19  REAOt IN.66)G.A.XM«ANG0.VI 
IF(G.LE.-I0.)G0T02 
REAOUN.206)W 
RRITC<IO.?7)W.P.TEMP 
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w«i  itno.  loiziSHAHt 
ikRITF  (lO.I013)XM 
VHUZ  * vi 

* XM/7000. 

*>jon  * *ngo  • .001 

Noncdo.  loi  1 1 VI 

pf,  PfAiM  Pi,?n  J)  VPLANt 
?fc  CONTlNOt 

hPIIt  VPLANt  < AI^GO 

VI  *IVI  ♦ VPLAVEI/3.2n0fi 
C»LL(.POUP  ( V I .VV.N.  1 > 

K»N- 1 » I 
00  J2*<  J * 

L«1*J-I 
C0(  J1*CCIU 
V(JI*VV(LI 
28  VtLDY  « 

V ( 1 I * VI 
T ( 1)  « 0. 

X ( n«fj. 

4fj(i  c A'.CiO 

r 1 » V ( 1 1 / 1 0 . 

r)vi»ovv/i'). 

lOVxOVl 

K*1 

C?i1.8h*A/I<,4./XK 


Vl)l  = 0 • 

OV*Dvv 

ZOl*!). 

Z X|l  • 

Kl-1-1 

K*«*10V 

KK«K-(IW 

30  YH»3,28l*Y ( l-l) 
ir(YH,Lf..O,  IGOTOM 

CALtrOOlOI YH,TEMP*PKE5>*0»V1S«VELA»0> 
r,0T032 

31  n«, 00237 

32  [)»J2, 1739*0 
TEHP«T£mP  - 459.69 
CaC2«0 

veAK*(V(l)  ♦ V(I-ll»*.5 
ANYjIxANG 
ANO«ANG- .0005 
ir(ANG.LT.-l.5IGOTO40 

221  CONTINUE 

0T«2.*(I./<v(n-Xl*UV»-I./<VU»”<Xl-l.I*DV))/C/(CD(l)*CD(I“l)) 

0T«  OT  » ?./C/(CDU)  » COn-IU*ALOG(COS(ANGl)/COS(ANG)  )/VBA« 

Av=.5* (C0(  I ) * ( V( 1 1 - XI*0V1**2*SIN(ANG)  ♦ CO ( l-l ) * ( V ( I ) - IX 

2l-1.0**OVI  **2“SIN(  ANGUl 
Till  all l-l I • OT 
VEUlYa  VOI  ♦ C*AV*0l/2. 

»NG2=ANG  , 

YOUTa  -VELOY  -0*T(l)  ♦ vm*SlN(ANOO)  - C*AV*DT/2. 

ASr,aAT*N(  YOOT/IVI  n-Xl*OV)/COSI  ANGI  ) 

IF(AHS(ANG2  - ANOI  .LE.  .OOI*APSIANG) )G0T0224 
0010221 

224  VOI*  VElOY  . C*AV*0T/2. 

222  , j , . I .|•0VI•C0S(ANGl)  - (Vll)  - Xl*DV  I *COS  ( ANGI  I / 

2C/I ICOI I I • I VI 1 1 - Xl*OV»l  ♦ ICOI l-l )*(V( I ) - (Xl-l.l*DVI  I I 

Y(i'i»Yn-i)  - !!5*(G*nm  ♦ tii-i)1)«dt  - veldy*dt  ♦dt»vi ii»sin( a 
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C»  (C0(  II»(V(  11“A1*0V)  ♦ CD  ( 1-1 ) • < V ( 1 1 - (Xl-1 . 1 *DV)  >*3.35 

ZD*  <<•  (ZfH  - KI»EAPC-CbA'»*OTl 

Z*  Z ♦ •'"01  ♦ 1301  - W»/CtJAM*(l.-EXP(-CBAR*DT>> 

Z0I*70 

73(11*/ 

ANOA  ( 1 I »AN(i 

end  i*c 

VELUA ( I > *VELUr 
IE  (Yd  I I I 1 • 13*13 
1?  1E(  «/d  I .LF.O.IGOIOl  1 
1F(  I. Of  .3nOIG')T011 

I n continue 

II  WUITM10.4A1 
1*1-1 
fjnli)j*l*l 

koITE. (10*‘t1I  Y(J»**(J)*T(J)«V<J)»CD(J)*77(J) 
in  CONTINUE 
kka*kKn  • 1 
r.orozi 

40  STOP 


ENII 

SUOnOUTINE  FOrtlO  (M«T*P*K.VK*VS*  IDPT) 

FOOIU  «‘‘S  PREVIOUSLY  AlMSl  

DIMENSION  H4(lRi.*l(lll*«(3dl>*W3dl)*TB(ll)*PB(ll)*WB(ll> 
niwFNSlON  A(T),H(J>tC(3).0(3) 

3Sn  format  ( si  HD  1 135a7S37<><»0  105160 170300700-335S69-S3561 3000  157609*/ 

1 6/Ml  38466 1 13(l<dOOO  130069- 159303-75931800  O3O63343A1450054 130006309* 

260M1 70034  0 754  141 34 16460  3507 156033960333 1399  76 1370051868831  16317*/ 
3S7H2376933899M847373O7O6238998BO5197977655O07aeOOO35155388O4* 

436h5O0708OOO131H1 1 39466398180*  / 

561H3 10841 1093901 06310094361406 1881556333333386 18075781 84535661 88. 
660950954 1 3 383036 10039 75961 13304  86356766490305763369801987391 000  * / 
734H759S1 1 1 741643300359357873739661814) 

351  FORMAT(ri.O*-3P6F3*0*5F3.0*  10PF7*0*  5PF 7 * 0 * 0P3E 1 * 0 * 9PF6*0*/ 

1 IOPF6*0*  4Rf 6*0«OP3F 1 ,0*  9PF6.0*  10PE7»0*  5PF7*0*0P3F 1 *0*  9PF6*0 

2 * 10PF6.0*  4RF5,0*0PF1 .0*  10PF6.0*  5PF6 * 0 * OPF 1 *0*  11PF6*0*  5PFb»0 

3 *0PEl*O*  11RE6.0*  5PF6.0*0PF 1 .0*  11PF6*0*  5PF6*0* OPF 1 . 0* 3P2F7 . 0 * / 

4 «Rf0.n*3P3E6*O*7PF4.O*3P3F6.O*8PF4*O*4P2F7*O*lOPF5*O* 

5 4R3F7*0* 1 1P76*0*3PF6.0*/5PF4,0*12PF5*0* 

6 3Mf6*0*7PF5*0*12PF4,0*  3PF6.0*ePF5*0*13PF4,0*3PF7*0*9PF4*0*15PF4*0 

7 * JPF7*0* 10PE5*0«15PF4,0*3PF7.0*  10PF5*0* 16PF4.0*4PF4*'0*0PF7*0* 

a 6PF0.O* 1 3RF6*0*2pF5.0*-4PF1.0*OP4F1.0*/6P2F6*0*-3P2F3*0*6P2F6*0* 

9 -4R2F2.0I  „ 

format  ( 1 X *F3.3*1  lF10.0*/23dK«3E30*10/)  *F8.4*F  1 0 . 0 *F  1 0 .6 . E20 . 8 * 

1 F0. 3*F1 0, 0*/l X*4< F5. 31/lX* 2F 10, 6* 2F 10*0*/ 1X*2E 10*6* 5A*2F 10.0) 

3001  format (1X.5HALT*  *E9.4*13M  BELOW  LIMIT) 

IF  (lTAPE.EU.4)fiO  TO  5 
1TAPE*4 
PFwlNO  ITAPE 
RP1TF(ITAPE*350) 

PE»IN()  ITAPE  . .. 

READ  (1  TAPE*  JSl)  (MHd)*l*l*12)*(Wl(l)*W2(l)*W3d>*l*l*ll)»(TH(I)* 
1 PH(1)*PB(1)*I*1*11)*CUN1 *CON3*CDN3*CON4*CON5*CON6* (A(1)*B(1)*C(I) 
3 *ud), 1*1*3) 

5 CONTINUE 

Hf,P*CONl»H/d.*(CONl"H/CON3)  ) 

IF(M0P.LT  .0.)  HGP*0. 

00  1003  M*l*ll 
IF (HGP-hh  <M) 1 1003*1004*1003 
1003  CONTINUE 

IF  ( (HGP-HH (13) ) .GT.O. » GO  TO  1053 
m*13 


1003  M*M-l 

1004  TM»T8(M)*d.»Rl<M)*<HGP-M(3(M)  ) ) 

IF  ((HGP-90000.1.GT.O.)  GO  TO  1006 
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T»TM 

r,0  TO  10  70 

100f>  IF  ( (HGP-la0000.»  .6T.0.J  GO  TO  1009 
U? 

60  TO  1007 
1009  lij 

1007  T = TM*(Am-B<l)*ATAN(<MGP-C(I)  1/0(1)  )) 

1070  IF  (M?(i4)  ) lOlltlOaOtlOll 
1011  TFMPsl .♦wl <M)»(mGP-H8(M) ) 

(M)  /TEHP«*W?(H) 

R=Wrt(M)  /TEMP*«<  l.*i<2<M)  , 

60  TO  1030 

1020  TEmP«EXP(-i»3(M1«<HGP-MH(M)  )) 

PsPH (M) *TEMP 
MrWH (H) *TEMP 

1030  IF  ( (HGP-rONO) .GT.O.)  GO  TO  1032 
VS«!C0N3*S0RT  (T>11 
VK*COM<**(T««l.S/(  <T*C0N5)*«)  ) 

1F(H.LT.0.)  go  to  1050 
MFIUWN 

1050  IF (lOPT.GT.O)  WRITE  < JOUTPT*2001 ) H 
RFTUR'^ 

105?  T*0. 

PiO  . 

Wsfl. 

103?  V5=0. 

VK*!). 

RETURN 

FNf) 

5UMR0UT  INE1NTER«<Z»ZZ»TT»T*J) 
niMENSlCNT (2001 *Z( 200) 

A*1  . 

1F(T(1) ,GT.T(2) )A*-1. 

3 1F(A*(TT-T(J) )) 112.100.113 
11?  ZZ*Z(J-1).((Z(J)-Z(J-1)  )/l,T(J)-T(J-l)))*lTT-T(J-l)  ) 
RETURN 
113  J*J»1 

60  ro  T 
100  ZZ«Z(J) 
return 

END 

SUHROUT INEGROUP(VHAH.V.N.K) 

D1MENS10NV(200) 

RaO 

CLOSEaA8S(V(H-VB*H) 

0051*1 ,N 

CL0SE*AM1N1  (Cl.0SE.ABS(V(l)-V9AR)  ) 
1F(CLOSE.LE.10.)GOT06 
5 CONTINUE 

5 1F(ABS(VBAR-V(I»1) ) . LT. CLOSE )KrK*l 
K*K»1 
RETURN 
ENO 

SUHROUT INFOROUP(VHAR.V.N.K) 

OIMENSIONV(200) 

KaO 

CLOSEsABS ( V ( 1 ) -VBAR) 

0051*1. N 

CLOSE* AM  INI (ClOSE.ABS(V(I)-V0AR)  ) 
1F(CLOSE.LE.IO.)GOTO(> 

5 CONTINUE 

6 lF(AflS(VBAR-V(I.l) ) .LT .CLOSE ) K*K.  1 
KsK*  1 

RETURN 

END 
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Graphs  1 to  140 
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fee  N«  wcs's  qm 


vn  M M9  «•  no  0 


26 


2Q  MM  RMMUN! TION 
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20  MT*I  fiMMUNiTION 

flIR  TO  GROUND  ftIRCRRFT  VELOCITT  790  KNOTS 

I Allltuoc  «)  HMCr 
II  iinr  *s  MNcr 

III  vriK'TT  vs  MMcr 

s«»i  }ioo  CMiiiS  ?tso.rfS  ruikc  viikitt 
^ USDS  ifso  &MIMS  Mso  rrs  NuiiLC  veiociTT 
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